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Abstract 

The problem of constructing global models describing isolated axially sym- 
metric rotating bodies in equilibrium is analyzed. It is claimed that, whenever the 
global spacetime is constructed by giving boundary data on the limiting surface 
of the body and integrating Einstein's equations both inside and ouside the body, 
the problem becomes overdetermined. Similarly, when the spacetime describing 
the interior of the body is explicitly given, the problem of finding the exterior 
vacuum solution becomes overdetermined. We discuss in detail the procedure to 
be followed in order to construct the exterior vacuum field created by a given but 
arbitrary distribution of matter. Finally, the uniqueness of the exterior vacuum 
gravitational field is proven by exploiting the harmonic map formulation of the 
vacuum equations and the boundary conditions prescribed from the matching. 



The description and analysis of isolated rotating objects in equilibrium is a fun- 
damental problem in general relativity. Unfortunately, despite important efforts in 
that direction, our present understanding of the question is still very incomplete. In 
particular, not a single model describing both the interior and the exterior of a self- 
gravitating, spatially compact, rotating object is explicitly known. There are, however, 
very interesting solutions which describe the exterior field of a physical object. We are 
referring to the solution for the rotating disc of dust recently found by Neugebauer 
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and Mainel This metric is vacuum everywhere and contains discontinuities in the 
first derivatives of the metric across an equatorial disc centered on the axis of sym- 
metry, which can be interpreted as a shell of dust. Despite its manifest interest and 
physical relevance, this solution does not contain an interior body in strict sense (in 
mathematical terms, the interior of the body is empty). In general, solving the vacuum 
field equations everywhere and allowing for discontinuities in the first derivatives of the 
metric will be adequate for describing the gravitational field generated by a shell of 
matter, or as a possibly good approximation of some very thin objects, but it cannot 
be used as an exact solution when the source object has a non-empty interior. In that 
case, the problem becomes even more difficult because the global model must be con- 
structed from two spacetimes, one describing the interior of the body (i.e. a solution 
of Einstein's field equations with matter) and another one describing the exterior field 
(i.e. an asymptotically flat solution of the vacuum equations). Several methods have 
been presented recently for attacking this kind of problem. One of them requires the 
resolution of Einstein's field equations both inside and outside the body starting from 
a common boundary. This approach is taken, for instance, in a recent paper where 
the existence, regularity and uniqueness of the solution has been investigated. More 
precisely, the authors prescribe Dirichlet boundary data on a two-sphere and then show 
that both the vacuum field equations and the interior problem (assuming rigid rota- 
tion) are solvable. The interest of these results is obvious, but it must be noted that 
they still do not prove the solvability of the global problem, not even for the rigidly 
rotating case. The reason is that it cannot be ensured that the two solutions thus 
constructed match appropriately across the boundary hypersurface. Given the type of 
boundary data used in 0, it follows that the global spacetime has indeed a continuous 
metric, but nothing can be said a priori about the first derivatives of the metric (they 
can certainly be discontinuous). Hence, in general, the global spacetime will have an 
undesired and uncontrolled surface layer of matter, which is an important weakness 
of the method. The same argument shows in general that, whenever two different 
spacetime pieces (interior and exterior) are used to construct the global manifold, the 
physical problem requires imposing the full set of matching conditions and not just a 
(mathematically suitable) subset of them. 

In order to solve the global problem one can alternatively adopt the point of view of 
finding the whole spacetime once and for all. The most favourable situation would be 
finding the equivalent of a Green function for the coupled system of equations describ- 
ing gravity and matter, so that the solution can be expressed directly in terms of the 
sources. Although this is probably the mathematically most well-motivated method, 
the technical difficulties are enormous and very scarce results are known at present. 
In particular, it seems very unlikely that an explicit model can be obtained using this 
method in the short, or even in the medium, term. It must be mentioned, however, 
that this approach has been followed in a very interesting paper by Heilig [^, where 
the existence of global solutions describing rigidly rotating stars is proven for configu- 
rations sufficiently close to a Newtonian solution. A second global method that can be 
employed requires finding smooth solutions of the differential system describing gravity 
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and matter, such that the matter fields tend to zero fast enough when we approach 
infinity (so that the solution can still be interpreted as an isolated body). Although 
this method is potentially capable of producing explicit solutions, its main drawback 
is that the rotating body is not truly finite in size because it has no limiting surface. 
Hence, the solutions have no vacuum exterior. Of course this does not mean that the 
solutions are necessarily unphysical, but for many purposes it is highly convenient to 
have models in which the rotating object has a clear limiting boundary. A second 
disadvantage shared by the two global methods described above is that the system of 
differential equations to be analysed must be closed, and this requires the imposition 
of conditions on the type of matter beforehand such as, for instance, rigid rotation and 
a barotropic equation of state. This means that a different analysis must be performed 
for each different type of matter. 

All in all, it follows that for any truly isolated and finite body, the full set of 
matching conditions between the interior metric and the exterior vacuum solution 
must be inevitably imposed in order to analyse the physical problem. Excluding the 
"Green function" approach, the construction of the global solution requires three steps, 
namely, obtaining the interior metric (either by solving the Einstein field equations for 
an adequate type of matter or by any other means), solving the exterior vacuum field 
equations and imposing the matching conditions between the two spacetimes. These 
three steps could be solved, in principle, in any order. In our opinion, however, the 
natural ordering is to provide first the interior metric and then finding the exterior met- 
ric which satisfies the matching conditions with the prescribed interior region. This 
point of view incorporates in a natural way the common idea that the exterior field 
is created by the interior distribution of matter. Moreover, the natural questions of 
existence and uniqueness of the gravitational field are only reasonable for the exterior 
problem, once the interior metric is known. The reason is that we cannot expect exis- 
tence, or even uniqueness, of the interior region for a given exterior gravitational field 
(it is well-known that in Newtonian gravity different distributions of matter produce 
the same exterior gravitational field; also, the spherically symmetric case proves that 
this still holds in General Relativity). Moreover, the inside-outside procedure has the 
important advantage that the resolution of the interior field equations, which are much 
more difficult to handle that the vacuum equations, can be performed by any method 
we wish, or even not done at all if we prescribe the interior metric by hand. In other 
words, there is no need to solve the interior equations in order to study the associated 
exterior problem. Obviously, for physically relevant problems the interior field equa- 
tions should be solved, but the important point here is that the interior problem is 
completely decoupled from the exterior one. 

The aim of this letter is, besides stressing and clarifying all the points above, to 
describe in detail the conditions one obtains by matching an interior (arbitrary but 
supposed to be known) stationary and axially symmetric metric with an unknown 
exterior vacuum metric, as well as to analyze and solve some of the problems and dif- 
ficulties that arise. In particular, we show that the junction conditions introduce two 
new essential parameters in the interior metric. To the best of our knowledge, this 
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fact was unknown hitherto. We also show that, for each value of these two parameters, 
the matching hypersurface, both from the exterior and the interior regions, is uniquely 
determined. Furthermore, the matching conditions fix the boundary conditions for the 
exterior problem. When written in terms of the Ernst potential the boundary value 
exterior problem turns out to be overdetermined and not unique. It is overdetermined 
because the interior geometry of the self-gravitating object provides more boundary 
conditions than necessary for solving the Ernst equations and it is not unique for, 
as we shall see, the boundary conditions fix the Ernst potential on the matching hy- 
persurface up to an arbitrary, non-trivial, additive constant. Since the Einstein field 
equations for stationary fields are of elliptic type, it follows that the physical problem 
is overdetermined and mathematically not well-posed. After computing and classify- 
ing the matching conditions, we finally prove the uniqueness of the exterior vacuum 
gravitational field generated by a given isolated, axially symmetric, rotating object 
in equilibrium. Although this result is quite natural and it has been often assumed 
implicitly in numerical calculations, no proof is available in the literature. The only 
results we are aware of are due to Lichnerowicz 0, but they only apply doubly locally, 
in the sense that they hold only in a neighbourhood of a local piece of the matching 
hypersurface. The method we use for proving the uniqueness rely on the harmonic 
map formulation of the exterior vacuum field equations. This provides a proof which 
is extremely simple and which is valid for any boundary data (by using more classi- 
cal methods of functional analysis one normally obtains theorems which only hold for 
boundary data satisfying certain inequalities). It must be emphasized here that none 
of the recent results on the uniqueness of solutions of the Ernst equations - see e.g. 
and references therein -, analyzes the problem we are considering in this letter. 

Let us then describe in detail the procedure discussed above for constructing global 
models for rotating objects. Let us assume we are given a stationary and axisymmetric 
spacetime (Vi, gj) describing the interior of the selfgravitating fluid. We make the usual 
assumption that this metric satisfies the circularity condition [0 (for fluids without 
energy flux, this is equivalent to the absence of convective motions). Hence, there 
exist coordinates {T, $, r, C,} in which the line-element reads 

ds] = -e^^ {dT + Bd<^f + [e^^ (dr^ + dC) + a'd^^] (1) 

where V, B, h and a are (known) functions of r and ( and 5$ is the axial Killing. 
The searched vacuum exterior manifold (Ve,^'^) is assumed truly stationary (free of 
ergoregions and/or Killing horizons) so that it can be globally described using Weyl 
coordinates 

dsl = -e^^ {dt + Ad(j)f + [e^^ + dz^) + P^dcf)^] (2) 

where U, A and k are functions of p and z only. The axial Killing vector of this 
spacetime is and the axis of symmetry is located at p = 0. The coordinates have 
been chosen so that dt is the unique timelike Killing vector which is unit at infinity (i.e. 
t mesures proper time at infinity) Then, the coordinate freedom in (^ consists only of 
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trivial constant shifts of t, (j) and z. Introducing the Ernst potential 0] £^ = e^^ + ifl, 
where Q is defined up to an additive constant by 

pr],, = -e^^A,, pfi, = e^^A,, (3) 

(a comma indicates partial derivative), the vacuum field equations read simply 

2AU + e~'^^{dn,dn) = 0, Ail - A{dn, dU) = , (4) 

where d is the exterior derivative, A is the laplacian of the 3-Euclidean metric 

dslj = dp^ + dz^ + 

and ( , ) denotes scalar product in this fiat space. The equation for the remaining 
function /c is a simple quadrature integrable once the solution of (|D is obtained. Thus, 
the system of differential equations to be solved constitutes a coupled, non-linear, 
quasi-linear elliptic system of p.d.e. where the principal part is the Laplace-Beltrami 
operator on the Euclidean 3-space. 

In order to find the boundary conditions of this problem, the matching of the 
known interior iVi^gi) with the undetermined exterior {VEigE) must be analyzed. The 
standard junction theory involves finding two isometric imbeddings of a 3-dimensional 
manifold (S, g^) into (V?, gi) and {Ve, gs) respectively such that the two second funda- 
mental forms in E inherited from both imbeddings coincide (see e.g. 0). It is natural 
to assume that S preserves both the stationarity and the axisymmetry and hence local 
coordinates {A, r, (/?} in S can be chosen so that the imbedding xe '■ S — > Ve reads 

XE-{t = r, (j) = (p, p = p(A), z = z{X)} , (5) 

that is, we adapt the coordinates in S so that dr and coincide with dt and on 
the exterior matching hypersurface xe{^)- The next task is the determination of the 
imbedding xi '■ S — > Vj. Obviously, the images of dr and d^p in xi (S) must be Killing 
vectors. Since the axial Killing vector has an intrinsic definition as the only one with 27r- 
periodic closed orbits and a symmetry axis, we must impose xi* (^(/^Is) ~ ^*lx7(s) (^i* 
is the push- forward of xi)- Regarding the image of dr, and given that the stationary 
Killing vector in {Vi,gi) does not have an intrinsic characterization, we must allow for 

XI* (drl^) = a (dr + bd^\^^j.-^ , (6) 

where a is a positive constant (to preserve the time orientation) and 6 is a constant 
restricted only to preserve the timelike character of the Killing. In order to deal with 
this freedom, we perform in {Vi,gi) the coordinate change 

$ = $' + abT\ T = aT' (7) 

which implies Ot' = a{dT + bd,i,) and 9$' = di^,, leaving the axial Killing invariant. 
The new metric potentials (defined so that (0) keeps the same form when all unprimed 
quantities are substituted by their primed counterparts) are a' = aa, h' — V = h — V , 
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Thus, the matching procedure has introduced two new essential parameters in the 
interior metric (even though they are irrelevant for the interior geometry by itself). 
In order to clarify this, let us show that matchings with different values of a and b 
are essentially inequivalent from the exterior point of view, by taking as an example 
an interior describing a non-convective fluid with velocity vector u = N {Ot + ^f(9$), 
where and uj are functions of r and C- For any matching satisfying (j^), w on the 
matching hypersurface as seen from the exterior becomes 



(we will identify T, = xi (S) = xe (S) from now on to simplify notation). The coor- 
dinates t and (j) in the asymptotically fiat exterior have an intrinsic meaning and (||) 
shows that the proper angular velocity of the fluid on S is a(w — fe), which depends 
essentially on a and b. Thus, different values of a and b will potentially give rise to dif- 
ferent exteriors (whenever they exist), because the state of motion of the interior body 
is distinct as seen from infinity. Consequently, we should not expect that the exterior 
gravitational field for different values of a and b coincides, and the physically relevant 
problem is the uniqueness of the exterior gravitational field when a and b are fixed. In 
other words, the exterior gravitational field will be fixed when the interior metric and 
the identification of the interior with the exterior through S are both prescribed. 

The new coordinate system (^ will be used from now on dropping primes every- 
where. The imbedding xi • ^ — ^ ^/ reads then 



The four unknowns p(A), z{X) in (^ and r(A), C(A) in (|) define the junction hypersur- 
faces and must be determined from the matching conditions. It can be shown [|I^ that 
all the matching equations can be reorganized into the following three sets of conditions 
(to prove the equivalence the vacuum field equations for the exterior have to be used): 
1 Conditions on the interior hypersurface. The functions r (A), ( (A) defining the inte- 
rior hypersurface are the solutions of the overdetermined system of ordinary differential 
equations 



6>o = {dT + Bd^) , 01 = ae-^d$, 6>2 = e^-^ dr, = e^-^dC, 



and the dot denotes derivative with respect to A. In many cases, these two differential 
equations are incompatible and no matching hypersurface exists, meaning that the 
given metric cannot be considered as a model for describing an isolated body with a 
vacuum exterior. There are also particular energy-momentum tensors for which the 





X/:{T = r,$ = (^,r = r(A),C = C(A)}. 



(9) 
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above equations do not provide any information at all (those with G22 = G23 = G^s = 0, 
including the important case of dust), so that the matching can be performed, in 
principle, at any timelike hypersurface preserving the symmetry. However, as is clear 
from ([TD|), the interior matching hypersurface is unique in the generic case. 
2 Exterior matching hypersurface. The exterior matching hypersurface p{X),z{X) is 
then uniquely determined by 



P(A) 



a 



z (A) = a^rC ~ ct,c'f' 



(the additive constant in z {X) is spurious given the shift freedom z — > z + const.). 
3 Boundary conditions for the exterior problem. Finally, the matching provides bound- 
ary conditions on the exterior metric potentials U and A at the matching hypersurface 



U\j,= V\^, 
U^z-U^p\^=VrC-V^r^, 



A 



B\ 



B^rC - B^^r 



(12) 



^From these equations we must get the appropriate boundary conditions for the exterior 
vacuum problem in terms of the Ernst potential. This can be done by noticing that 



n 



= Q^pp + Q^J 
Q^pZ — ^l^zP 



— {-A^,p + A^pz] 



P 



(A ,i + A^pp) 



a 



rv \ 

{B,.r + E^C 



where the righthand sides are known. Thus, the matching conditions fix the normal 
derivative of on S uniquely, but they only fix O on S up to an arbitrary additive 
constant. As we shall see, this fact will add some subtleties in the uniqueness proof 
below. 

We are left with Ernst equations (^) subject to boundary conditions both at the 
fixed (and known) compact axially symmetric surface given by (|ll]) and at infinity. 
The conditions at infinity ensure that the spacetime is asymptotically flat (which is 
the mathematical translation of the isolation of the body) and are given by 



U = -niR 



-1 



0{R- 



-2zJR-^ + 0{R~ 



(13) 



where R = -y/p^ + 2;^ and m and J are the total mass and angular momentum of the 
source, respectively. It is now evident that the boundary conditions on S constitute a 
Cauchy problem for the elliptic system (|^), so that it is mathematically not well-posed. 
Thus, the physical problem is overdetermined and we should not expect existence of 
solutions in all situations. This is the main problem for proving the existence of the 
exterior solution given the interior, and it is a matter of further investigation to flnd 
which restrictions on {Vj,gj) will allow to overcome it. 

Nevertheless, we can now address the proof of the uniqueness of the exterior fleld, 
which is one of the purposes of this letter, by using the theory of harmonic maps. 



7 



Roughly speaking, harmonic maps are apphcations between Riemannian manifolds 
that extremize the so-called energy or Dirichlet functional. More precisely, given an 
n-dimensional Riemannian manifold (possibly with boundary) (M , ^) , a A;- dimensional 
Riemannian manifold (A^, 7), and a map \I' : {M,g) — > (A^, 7), we can construct 
the energy functional = /^^ ^ [tvM'^^i'^)] f?, where 77 is the volume form of M 

and trjvf\l/*(7) is the trace in M of the pull-back of the metric tensor 7 of the target 
manifold A^. If we introduce local coordinates {x^} and in M and A^ respectively, 
then becomes 

where /x, run from 1 to n and a, h run from 1 to k. Then, any map with finite 
E is called harmonic iff it is a critical point of E with respect to variations vanishing 
on the boundary of M (whenever dM 7^ 0). Thus, any harmonic map \E' satisfies the 
corresponding Euler-Lagrange equations, which read (in local coordinates for A^) 

A^" + r^;,(^(x))(£^^", c^^^) = 0, (14) 

where F^^ are the Christoffel symbols in A^. It is well-known that Ernst equations (^) 
can be written in this Euler-Lagrange form if the domain manifold M is chosen to 
be the part of the Euclidean 3-space with extends from the two-dimensional compact 
boundary S = {p(A), 2;(A), 0} to infinity and the target manifold A^ to be the two- 
dimensional manifold with metric 

dsl = AdU^ + e-^^dn^ (15) 

(Notice that S is in fact the image of a new imbedding of the abstract S, in this case 
into M. Thus, even though there is a one-to-one correspondence between S and S, 
they should not be confused as their corresponding first fundamental forms do not 
coincide). The solutions of (^ we are interested in must must have the asymptotic 
behaviour (|13D which ensures that the corresponding energy functional is finite. 

The proof we will give for the uniqueness of the exterior gravitational field consists 
in two steps. First, we will show that two solutions with the same Dirichlet boundary 
conditions on S must coincide. To that aim, let us assume that there are two harmonic 
maps \l/o and from M to A^ satisfying Uo\j. = Ui\j. , fiols = To proof 

that they coincide everywhere, we will make use of a most remarkable identity due to 
Bunting |Tl|]. The idea is to construct a smooth homotopic application 



^ : M X [0, 1] — > N 

{x,u) — > -^{x,u) (16) 

satisfying \l/(a;,0) = \I'o, \E'(x, 1) = \I'i everywhere on M and such that jxo{u) = 
\E'(xo,m) is an affinely parametrized geodesic in A^ for every xq G M. In our case, A^ 
is a simply connected geodesically complete maximally symmetric plane with constant 
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negative curvature so that a well-known theorem (see e.g. |T2|) asserts that any two 
points can be joined by a unique geodesic. Therefore the homotopy (|TB]) can always be 
constructed. A remarkably simple calculation (see also |jl3[) leads to Bunting's identity 



iv^V^s^ = J^' du [V'^s^V.s' - i?,,Hs'V^^^s'=V^^'} (17) 

where Rijki is the Riemann tensor in A^, s*(x) is the tangent vector along the geodesic 
7a; and s^(x) is its norm (being geodesic and affinely parametrized does not depend 
on u). For with non-positive curvature, the righthand side of this identity is non- 
positive definite and it only vanishes when s* is covariantly constant. The lefthand side 
in (plTl) is a total divergence and its integral in M can be transformed into a surface 
integral on the boundary S and the integral "at infinity" . The coincidence of and 
on E implies s|g = and consequently the surface integral on S vanishes. Regarding 
the integral "at infinity" , the asymptotic behaviour of \l/o and is of the form (|13D 

= (-rriDR'^ + O (r-^) , -2Jd zR'^ + O (R~^)) , D = 0,1 



(we cannot assume a priori the same mass and angular momentum for both solutions). 
The geodesies of A^ can be explicitly integrated so that the evaluation of the geodesic 
distance between \&o(a^) and \l/i(x) at points x lying on a spherical surface R = Rq 
proves that s falls off much faster than R^^ when Rq tends to infinity. Hence, the 
integral on the surface "at infinity" also vanishes and we can conclude that is a 
constant vector field in M. Since s* vanishes on S, s* must be zero everywhere and 
\l/o(a^) = ^^lix) Vx G M as we wanted to prove. 

In the second part of the proof we must take into account that, as shown above, 
the matching conditions fix the value of f2 on S (and thereby on S) up to an arbitrary 
additive constant which is relevant , because the freedom in Q has already been used to 
set the asymptotic behaviour ([TB|). Therefore, we still need to show that this arbitrary 



additive constant does not introduce a one-parameter family of solutions. To that end, 
the full set of boundary conditions {\L2^ on U and Q will now help. Let us assume there 
exist two solutions = ([/+,r2+), and \E'_ = (f/„,f2_) satisfying 

u+\^ = u.\^, v^[/+|g = V^[/_|5, n+\^ = n.\^ + c, v ^n+\^ = v ^n_\^ (is) 

where C is an arbitrary constant. A direct calculation shows that, given any solution 
\1/ = {U, Q) of the Ernst equations (^), the three-parameter one-form on M 



C3 



(in 



W (^; ci, C2, C3) = 2 (C2 + 2c3fi) dU + [e"'^ (ci + cafi + 
has vanishing divergence, (i.e. V^W^ (\E') = 0) and, consequently, 

l'Wf,dSf'= j ly^d^'^ = 87rc2m(^) (19) 
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where 6*00 stands for the surface "at infinity" and the last equahty follows immediately 
from the asymptotic flatness conditions (0). On the other hand, relations (|TBp imply 

W{^+;Ci,C2,C3)k = VF($_;Ci+C2C + C3C^C2 + 2c3C,C3)|- . 

The combination of this expression with the identity (p!9D produces 

87rc2m(*+) = Stt (c2 + 2csC) m(^_) 

and, given that C2 and C3 are arbitrary constants, it necessarily follows 

m(^+)=m(^_), m(^_)C = 0. 

Since for physically well-behaved solutions the total mass cannot vanish, C must be zero 
and the uniqueness of the exterior solution generated by a given interior distribution 
of matter in stationary and axially symmetric rotation is completely proven. 
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